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Abstract 



In a recent paper entitled "High temperature expansion applied to 
fermions near Feshbach resonance" , (Phys. Rev. Lett. 92 160404 (2004)), Ho 
and Mueller have demonstrated a remarkable similarity between its high and 
low temerature properties at resonance. The quantum second virial coefficient 
plays a crucial role in their analysis, and has a universal value at resonance. 
In this paper, we explore the connection between the quantum and classical 
second virial coefficients, and show that near a Feshbach resonance an exact 
mapping from the quantum to classical form is possible. This gives rise to a 
scale-independent inverse square effective potential for the classical virial co- 
efficient. It is suggested that this may be tested by measuring the isothermal 
compressibility of the gas on the repulsive side of the resonance.. 
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Great advances have been made in the study of ultracold trapped atoms in recent times. 
The interaction between the atoms is normally weak in a dilute gas, but may be enhanced 
drastically by making use of the so-called Feshbach resonance [1] . This may be achieved by 
applying a magnetic field to tune the energy of the Zeeman level of the trapped atoms to one 
of the molecular resonances. At a zero-energy resonance, the scattering phase shift is n/2, 
and the cross section reaches the unitary limit. It is argued that the properties of a gas at 
Feshbach resonance have a universal character [2], irrespective of the type of gas in the trap. 
In this connection, Ho and Mueller [3] recently considered the interaction energy of a two- 
component Fermi gas near a Feshbach resonance. They compared the available experimental 
data of the ultracold gas in the vicinity of a Feshbach resonance with the predictions in the 
Boltzmann regime using the quantum second virial coefficient. Surprisingly, it was found 
that the high temperature expression, when used at micro-Kelvin temperatures, explains the 
data reasonably well. In the present work, we examine the use of the quantal second virial 
coefficient in more detail, and its connection to its classical counterpart at the Feshbach 
resonance. We find that the quantal expression may be mapped onto the classical formula 
exactly in the vicinity of the Feshbach resonance using an inverse square interaction in the 
configuration space. 

To appreciate this in more detail, note that the virial expansion of a (one-component) 
gas (classical as well as quantal) for the pressure P at a temperature T is given by [4] 

— = l + a 2 (nA 3 ) + a 3 (nA 3 ) 2 + .... , (1) 
nr 

where 03... etc are the dimensionless second, third. .virial coefficients, r = /c^T, n = N/V 

is the number density of particles, and A = \J ^nfi 2 / vtlt is the thermal de Broglie wavelength. 
The free energy F = E — TS may be easily obtained by integrating the above P with respect 
to the volume V (since P = —(dF/dV) T ), and hence also the energy E = — t 2 [<9(F/t)/<9t]. 
After subtracting out the energy of the perfect gas part (classical or quantal), one obtains 
the virial series for the interaction energy : 
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(2) 

This expression is valid for a classical as well as a quantum gas, provided the appropriate 
classical or quantal virial coefficients are used. Whereas for a classical gas the virial coef- 
ficients may be expressed as integrals involving the interaction potential, for the quantum 
problem a solution of the j-body problem is needed to obtain aj. Henceforth we shall denote 
the classical virial coefficients by Aj to differentiate from their quantum counterparts, to 
be still denoted as aj. For an interesting example where classical and quantum results are 
very different, consider a hard sphere classical Boltzmann gas. For a hard-sphere diameter 
f'a A2 = ^f(5f ) 3 > an d the higher order A/s may be expressed as Aj = k(A2Y~ 1 , where k is 
a constant. Substituting this in Eq.(2), we see that the interaction energy of the classical 
hard sphere gas vanishes identically in each power of (nA 3 ). This, of course, is not the case 
for a bosonic or fermionic quantum gas at low temperatures. 

We note that for a classical gas, (nA 3 ) << 1, and the series (2) may be terminated after 
the first order term in most cases: 
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The classical second virial coefficient A 2 is given by [4] 
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where (3 = 1/r, and V(r) is the interatomic two-body potential. The quantum second 
virial coefficient requires a knowledge of the bound-state and the continuum spectra. For a 
one-component quantum gas, b 2 = — ct2 is given by [4] 
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In the above, E n /s are the two-body bound states (if any), and 5i is the scattering phase 
shift in the /th partial wave. In the quantum case, for a Fermi gas of atoms interacting 
with a delta-function pseudo-potential, the spin-up atoms interact only with the spin-down 
ones, and one may consider a two-component (up, down) gas with = = n/2. In this 
situation, n 2 in (3) should be replaced by n 2 /4. Often, the multiplicative prefactor in Eq.(6) 
is omitted, being transferred instead to the RHS of Eq.(3). The modified equation for a 
two-component quantum gas then reads 
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Ho and Mueller [3] used the above formula for the quantum virial coefficient in the high 
temperature expression (7) to explain the low-temperature data. Two questions immedi- 
ately arise regarding this treatment : 

(a) For the quantum gas, since (nA 3 ) is larger than, or of order unity, higher order terms in 
Eq.(2) may be important. Presumably, near the Feshbach resonance, it is still the two-body 
term that dominates the interaction energy. 

(b) One may also argue that while using the high-temperature expression (7) for the inter- 
action energy, one should take the classical Eq.(4) for the second virial coefficient coefficient, 
instead of Eq.(8). 

In the following, we attempt to answer the second question in some detail. Note that at 
low temperatures (large f3), only small values of k contribute dominantly in the integrand 
in Eq. (6), and therefore only the / = partial wave in the sum is important. Therefore, 
rather than the total classical A2 given by Eq. (4), we should consider its I — classical 
component. Making the decomposition 
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it may be verified easily (by replacing the sum over / by an integral), that 
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One may enquire about the relation between the "semiclassical" expression (10) and the 
quantam counterpart — b-p ■ This question was actually answered nearly seventy years back 
by Kahn [5] in his doctoral dissertation. We shall give here a new compact derivation of his 
result, and for simplicity consider a potential with no bound states present. Using (8), and 
the relation E = 0_ we have 
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With this definition of the virial coefficient, at a zero energy resonance, = 1/2 or —1/2, 
according as the scattering length is negative (attractive interaction), or positive (repulsive 
interaction). We also note that b^iP) is just the Laplace transform of the derivative of the 
phase shift with rspect to the energy. In the lowest order WKB approximation, the phase 
shift is given by [6] 
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In the above, the effective potential Vi(r) is defined as 
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In Eq.(12), ri and r are the classical turning points where the respective integrands go to 
zero. The derivative of the phase shift may be written as 
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in terms of the unit step function Q(x). By noting that the Laplace transform of ^/§=^ is 
y^exp(— (3/ji), we immediately obtain the desired result 
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The multiplicative factor of 4^ may be understood by noting that b$ as defined by Eq.(8), is 
the quantum canonical two-body partition function for the relative motion in the I th partial 
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wave, without the noninteracting continuum part. Noting that the radial motion in a given 
partial wave is one dimensional, the classical limit of is 
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If the lowest order WKB approximation (14) for a given potential were exact, the classical 
and quantum (partial wave) virial coefficients would be related by = —A^/V^- But 
in general, specially near the Feshbach resonance, the WKB approximation for the phase 
shift is not accurate. Moreover, the universality property of the quantum b^ = ±1/2 at 
the resonance is not reflected by its classical or WKB counterpart. The latter invariably 
depends on temperature as well as the parameters of the potential (even at resonance). As 
an illustration, consider an attractive square-well potential of depth — V , and range R for 
the two-body interaction. When the dimensionless "strength factor" s = ^^^P-R 2 is unity, 

there is an s-wave zero-energy bound state, with b^ — |. If, on the other hand, we were to 

use the WKB expression (15), we would obtain \b^]wKB = "^( ex P(/^K)) — 1)- Unlike the 
exact quantum result, this depends both on temperature and the potential parameters, so 
that universality is lost. 

The only exception to this rule is for a scale-independent inverse square potential of the 
type 



where a is a dimensionless (real) parameter. It is known that the WKB method gives the 
exact phase shift [7] for such a potential, and it is easy to check that the classical virial 
coefficient, evaluated using Eq.(lO), is 



Using Eq.(17) we see that for a 2 = 1, we get b 2 = —1/2, as is appropriate for a repulsive 
potential at a zero-energy resonance. 

We now show that the classical form of the second virial coefficient may also be recovered 
from the quantum expression through an exact mapping near the zero-energy resonance. To 
this end, we start with Eq.(ll) for the I = partial wave, 



Near zero-energy (k — > 0), we may use the shape- independent effective range expansion, and 
for I — it is given by 
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In the above, a is the scattering length that approaches ±00 near a resonance, and r is the 
effective range. The virial coefficient near the resonance is governed by a alone, and we may 
ignore the second term. In that case, using (20), we get 



This may be rewritten as 
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By making the transformation k = a/x, we get 
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We have already noted, from Eq.(19), that the second term on the RHS may be reproduced 
classically by an inverse square potential with a — 1 for the repulsive case. The first term 
on the RHS, on performing the integral, is independent of a, and gives a correction of order 
h to the classical value in the neighbourhood of the resonance. We may thus write the above 
equation in the form 



6 (o) 



sgn(a) 



(-*) 


\ V2X 


\ / res 


\a\ir 



(26) 



where ^l^ '' is obtained from Eq.(lO) with the inverse square potential (18) of strength a — 1. 

We have shown above that as the Feshbach resonance is approached from the repulsive 
side, the two-body physics may be simulated by a universal inverse square potential of unit 
strength. This does not appear to hold on the attractive side. The inverse square potential 
in one space dimension is known as the Calogero-Sutherland potential [8], and has been 
thoroughly studied, including its two-point correlation function [9]. The integral of the 
latter is related to the isothermal compressibility, and can be computed numerically [10]. Its 
experimental measurement near the Feshbach resonance may be of interest, specially in view 
of the findings in this paper. The other intriguing point is the so-called Efimov effect [11]. 
When the two-body scattering length approaches infinity from the attractive side, a very 
large number of three-body bound states appear. More over, Efimov has shown that the 
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universal effective three-body interaction in this case is — 2 ^ R2 , w±±c±g u — ^v' 12 '23 
r| 1 )/3, and Sq ~ 1. We find, from the two-body analysis, that an inverse square two-body 
effective potential also arises at Feshbach resonance, albeit on the repulsive side. More 
research in this direction is planned. 
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